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A relativist ically invariant quantum theory first advanced by Bakamjian and Thomas has proven 
very useful in modeling few-body systems. For three particles or more, this approach is known 
formally to fail the constraint of cluster separability, whereby symmetries and conservation laws 
that hold for a system of particles also hold for isolated subsystems. Cluster separability can be 
restored by means of a recursive construction using unitary transformations, but implementation is 
difficult in practice, and the quantitative extent to which the Bakamjian- Thomas approach violates 
cluster separability has never been tested. This paper provides such a test by means of a model of a 
scalar probe in a three-particle system for which (1) it is simple enough that there is a straightforward 
solution that satisfies Poincare invariance and cluster separability, and (2) one can also apply the 
Bakamjian- Thomas approach. The difference between these calculations provides a measure of the 
size of the corrections from the Sokolov construction that are needed to restore cluster properties. 
Our estimates suggest that, in models based on nucleon degrees of freedom, the corrections that 
restore cluster properties are too small to effect calculations of observables. 

PACS numbers: 21.45+v 

I. INTRODUCTION 

There are two distinct requirements for describing quantum mechanical systems of particles under the requirements 
of special relativity. The first requirement is Poincare invariance: probabilities, expectation values and ensemble 
averages for equivalent experiments performed in different inertial frames are identical. A necessary and sufficient 
condition for a quantum theory to be Poincare invariant is that the dynamics is described by a unitary representation 
of the Poincare group [1]. The second requirement is cluster separability: isolated subsystems must have the same 
observable properties as they would in a framework in which the other "spectator" particles are absent entirely. 
This requirement justifies tests of special relativity on isolated subsystems. It applies both to systems of particles 
interacting among themselves (e.g. via the strong interaction) and to the current operators that allow them to interact 
with external fields. 

There is more than one way to implement the requirements of Poincare invariance in a quantum system, and cluster 
separability is not necessarily an automatic consequence of the implementation. In local quantum field theory, Poincare 
invariance and cluster separability are satisfied formally as a consequence of the covariance, spectral properties and the 
local commutation relations of the field operators. For systems amenable to perturbation methods, these properties 
can still hold. For non-perturbative systems (such as those based on the strong interaction) controlled approximations 
are replaced by a truncation or an expansion scheme, which replaces an infinite hierarchy of coupled non-linear field 
equations with a more tractable finite subset, and which may be valid only in restricted parameter domains. These 
schemes may still exhibit Lorentz covariance, but cluster separability does not automatically follow, and must be 
demonstrated on a case-by-case basis that depends on the scheme. 

An alternative implementation of Poincare invariance in a few-body system involves direct construction of the 
dynamical representations of the Poincare group in the presence interactions with a finite number of particle (rather 
than field) degrees of freedom. Such a construction was originally studied by Dirac [2], and Bakamjian and Thomas [3] 
provided an explicit construction of the generators for a system of two interacting particles. 

Many realistic, Poincare invariant, quantum mechanical models of strongly interacting few-body systems are based 
on a generalization of the Bakamjian-Thomas construction [4]. This construction can be used to construct models of 
systems of arbitrary numbers of particles and systems that do not conserve particle number. A representative sample 
of few-body applications of Poincare invariant quantum mechanics based on the Bakamjian-Thomas construction in- 
cludes relativistic constituent quark models [5-9], relativistic few-nucleon models [10-14], relativistic models involving 
electromagnetic probes [15-17], and relativistic models with particle production [18, 19]. 

The virtue of the Bakamjian-Thomas construction is that it provides a means for constructing Poincare generators 
for systems of interacting particles, and where the problems have a finite number of degrees of freedom and are in 
principle solvable. The limitation of this construction is that for systems of more than two particles, the dynamical 
representation of the Poincare group does not become a tensor product on states representing asymptotically separated 
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subsystems. While this limitation does not lead to observable consequences in the three-body S-matrix or bound-state 
observables, there are observable consequences when the three-body system is embedded in the four-particle Hilbert 
space, as it is in four-body problems. 

The most ambitious applications of the Bakamjian- Thomas construction are relativistic Faddeev calculations [10- 
14] of three-nucleon scattering, which have been performed using realistic nucleon-nucleon forces and three-nucleon 
forces. The next step in developing a relativistic few-body theory would be to model four-nucleon systems or electron 
scattering from a three-nucleon system. In both of these applications a three-nucleon subsystem is embedded in a 
four-particle Hilbert space, and, for the first time, there is the possibility of observable consequences of violations of 
cluster properties. 

The restoration of cluster properties can be achieved through a recursive construction due to Sokolov [20, 21]. 
The Sokolov construction generates a dynamical representation of the Poincare group satisfying cluster properties 
using of a hierarchy of unitary transformations, each of which preserves the S matrix while transforming tensor 
products of subsystem representations of the Poincare group to representations where the interactions can be added 
in a manner that preserves the underlying Poincare symmetry. The Bakamjian-Thomas representation is retrieved 
by setting the relevant unitary transformations to the identity operator. The Sokolov construction is non-trivial and 
the unitary transformations depend upon the interactions. An isolated three-body system is a special case where the 
cluster property can be achieved by means of a single overall unitary transformation that preserves the S matrix. For 
systems of more than three particles there are observable differences between the S matrices in the two representations. 
The interested reader can find a complete discussion in [21]. 

The construction of the Sokolov hierarchy of unitary transformations is sufficiently complicated that the technique 
has never been used in realistic calculations. Before undertaking computationally intensive four-body calculations, 
one would like to know the importance of the corrections required by cluster properties. While Ref. [21] argued that 
these corrections should be small in nuclear physics applications, this was never quantified in any model calculations. 
An additional investigation is needed to determine if these corrections can be ignored, can be treated perturbatively, 
or must be treated exactly. To address this question we construct a model involving matrix elements of a scalar probe 
for a system of three particles. The model is simple enough that both Poincare invariance and cluster separability 
are easily satisfied; this unusual pair of features is due to the simplicity of the test model. At the same time the 
model is sufficiently rich that it permits an equivalent Bakamjian-Thomas treatment of the three-body system that 
illustrates the quantitative impact of the breakdown of cluster separability in this four-body problem. The conclusion 
of our preliminary analysis suggests that the corrections required by cluster properties are too small to be observable 
in nuclear physics applications. The corrections are more important for models based on sub-nucleon degrees of 
freedom. 



In a few-body quantum mechanical model the requirement of cluster separability means that the Poincare generators 
{H, P,J,K} that generate time translations, translations, rotations, and rotationless Lorentz transformations and 
current operators I ti {x) have cluster expansions of the form 



II. CLUSTER PROPERTIES 




(2.1) 
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where the terms Xjj, Xijk, ■ ■ ■ are short-range 2, 3 • • • -body operators (which only depend on the degrees of freedom 
of the particles in a cluster and vanish when one particle in a cluster is separated from the remaining sub-cluster), 
and the sums of the generators and currents over any subset of particles satisfies the Poincare commutation relations, 
current covariance, and current conservation. 

These constraints are realized in the Sokolov construction, but they are not generally satisfied by the Bakamjian- 
Thomas construction. The Sokolov construction generates all of the many-body operators that appear in these cluster 
expansions as functions of input two-body interactions. These interactions are necessary to preserve the Poincare 
commutation relations in the presence of interactions. 

III. TEST MODEL 

Poincare invariance and the Dirac forms of dynamics for Bakamian- Thomas (BT) constructions are discussed 
extensively in Ref. [4]. The framework provided here makes use of that discussion. Conceptually, the Bakamjian- 
Thomas construction defines the dynamics in the rest frame of the non-interacting system and then uses simultaneous 
cigenstates of the mass and spin (Casimir operators for the Poincare group) to construct unitary representations of 
the Poincare group. We approach the model in a heuristic way in order to illustrate the issue of cluster separability 
combined with Poincare invariance with a minimum of formal development. 

A simple four-body model where we can observe the breakdown of cluster separability consists of a system of three 
spinlcss particles interacting with an external probe that represents the presence of a fourth particle. The relevant 
dynamical quantities are matrix elements of a one-body scalar density between eigenstates of the three-body four- 
momentum operator. We will see that, even for this simple model, there are differences in matrix elements of the 
scalar density calculated using three-body eigenstates of a BT four-momentum operator compared to matrix elements 
calculated using three-body eigenstates of a four-momentum operator that clusters. Our simplifying assumptions 
show that the violation of cluster properties has nothing to do with particle spins, or properties of more realistic 
probes like four-current densities. 

The scalar density operator, j(x), is assumed to be a one-body operator. The three-particle system consists of a 
bound pair (which we label the (12) subsystem) and a third particle (3 which we call the struck particle) that does 
not interact with particles 1 or 2. Only particle 3 interacts with the scalar density j(x). Therefore, the (12) subsystem 
acts as a spectator with respect to the action of the probe. This can be illustrated by disconnected graph shown in 
Fig. 1. 

Ps P12 



P3 Pl2 

FIG. 1: Graph of one-body density interacting with a particle plus a bound two-body spectator. 

In the following sections we will calculate matrix elements of the density operator in this simple system using 
the Bakamjian-Thomas unitary representation of the Poincare group, and compare them with a result that exhibits 
manifest cluster separability. The latter result is possible only because of the simplicity of the model. 

A. Tensor-Product (TP) Representation 

For the simple system introduced above, the unitary representation of the Poincare group is a tensor product of 
a two-body representation on a two-particle Hilbert space and a one-particle representation acting on a one-body 
Hilbert space. In this section we consider a three-body eigenstate that is a tensor product of bound (12) eigenstate 
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and a one-particle eigenstate 

|A,m;pi2,p 3 ) TP = |A;pi 2 ) ® |m,p 3 ). (3.1) 

where A labels the mass eigenvalue M\ of the bound state of the interacting two-body mass operator, M 12 , P12 is its 
three-momentum, and m and p 3 are the mass and momentum of particle 3. Since the scalar density j(x) acts only in 
the space of particle 3, the matrix element has the form 

Tp(A,TO;p' 12 ,p^|j(0)|A,m;pi2,p 3 ) T p = S(p' 12 - pi 2 )(m; p 3 |j(0)|m; p 3 ), (3.2) 

where 

\ju m (p 3 )(m,p 3 \j(0)\m,p 3 }y/u m (p 3 ) = mG(q 2 ); q = p' 3 -p 3 , (3.3) 

and q is the four- momentum transferred to particle 3 and w m (p 3 ) = \/m 2 + p§. The square roots in Eq. (3.3) ensure 
the Lorentz invariance of the form factor G(q 2 ) when the states of particle 3 have a delta-function normalization: 

(to, p^m, p 3 ) = 5(p 3 - p 3 ). (3.4) 

If we specify initial momenta, p 3 and pi 2 , and the momentum transfer, q, and integrate the matrix element over 
the final spectator momentum, p' 12 , the integral collapses due to the spectator momentum delta function in Eq. (3.2), 
and we have 

Ftp ■= j dp' 12 ^Juj m (p' 3 ) TP (X, m; p' 12 , p' 3 |j(0)|A, to; p i2 , p 3 ) TP V w m(P3) = mG(q 2 ). (3.5) 

Equation (3.5) is the Tensor-Product (TP) result. It has no dependence upon the momentum pi 2 of the bound-state 
spectator, as expected from the physical requirement of cluster separability, and does not depend upon the specific 
values of p 3 and p' 3 , as long as p' 3 — p 3 = q. The result is just the invariant form factor of particle 3 multiplied by its 
mass. 

The representation for a three-particle system becomes a tensor product of a one-body and two-body representation 
of the Poincare group whenever the interactions involving one of the particles is set to zero. This implies that 
the generators have the form (2.5). Ideally one would simply perform all calculations using the TP representation. 
However, for examples beyond the simple model presented here, it is a highly non-trivial task to obtain a tensor-product 
representation that preserves Poincare invariance. The only known method for systems with a finite number of degrees 
of freedom is to obtain first the Bakamjian-Thomas representation (which has interacting generators that satisfy the 
Poincare commutation relations), and then to restore cluster separability by a recursive construction involving a 
hierarchy Sokolov's unitary transformations. In the special case of the three-body system the two representations are 
related by a single unitary transformation, A: 

\iI))tp = A|V)bt. (3.6) 

The operator A is interaction dependent and contains the corrections that restore cluster properties to the BT 
representation. The three-body TP and BT representations related by A have the same 3-body S matrix, but off- 
shell differences lead to different predictions in the four-body system. The operator A converts a BT representation to 
a TP representation satisfying the constraints imposed by the Poincare commutation relations. However, for realistic 
three-body systems, A is a complicated interaction-dependent operator. 
The size of the difference 

tp(^'\j(0Wtp -bt W\j(0Mbt, (3.7) 

where the BT and TP states are related by Eq. (3.6), provides a measure of the importance of the corrections needed 
to restore cluster separability. 

While in general the Sokolov operators A need to be computed, the test model presented here is simple enough 
that the TP and BT states related by A can be computed directly, and then used in (3.7) to estimate the size of the 
corrections needed to restore cluster separability. 
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B. BT Representations 

In the following sections, we evaluate matrix elements of the scalar density (3.2) using BT representations of the 
2+1-particle system corresponding to different forms of dynamics [2] and compare these results to the tensor product 
result. These matrix elements will each depend upon the momentum dependence and mass eigenvalue of the two- 
body bound-state wave function, in violation of cluster separability, and our goal is to examine the magnitude of that 
violation. In all cases, we specify the four momentum q transferred to the struck particle and bound-state momentum 
P12 for a given calculation and then vary these momenta for sensitivity tests. Dependence of this matrix clement on 
P12 indicates a failure of cluster properties. 

In addition to the external momenta p 3 and pi 2 , the calculations depend upon the momenta of the constituent 
particles, which we label bi , b 2 and b3. The constituent momenta are not experimentally accessible, and their relation 
to the external momenta depends in turn upon the form of dynamics, as is discussed in detail below. For purposes of 
comparison we vary the external momenta keeping h' 3 — b 3 = q. This requires different kinematic conditions for each 
form of dynamics. 



C. BT Representation: Instant Form 



We now evaluate the scalar density matrix element in Eq. (3.2) using an instant-form BT model. Instant-form 
models have no interactions in the generators of rotations and space translations. 

We provide first an observation about frames that will be relevant to each of the BT representations. Calculations 
of current-matrix observables typically use a Breit frame (in which the energy transfer is zero), a lab frame (initial 
target three-momentum is zero), or an anti-lab frame (final target momentum is zero). There is no frame in which 
both initial and final target three-momenta can be zero, so the calculation will require Lorentz transformations relating 
frames with different target momentum. These are provided explicitly where needed below. 

First, we change variables, replacing the bound state and particle 3 momenta by the total momentum, P, of 
the system and the momentum of particle 3 transformed to the frame in which P = by a rotationless Lorentz 
transformation A _1 (P/M): 



where 



and 



(Pl2,P3) (P,P), 



P:=Pi2 + Ps; p:= A-^P/MJpa; p = p 3 + $_(p 3 , P, M)P, 



(3.8) 



(3.9) 



$_(p 3 ,P,M) 
M 

Wm(P3) 



1 

M 



P-P3 



Wm(P3) 



E + M 



With this variable change the relation between the (12) bound state and particle 3 in these bases is: 

0(P,p) 



|A,m;pi2,p 3 ) 



<9(P12,P3) 



|A,m;P,p), 



(3.10) 



(3.11) 



where the suffix A in the Jacobian indicates that the interacting two-body mass eigenvalue M\ was used in the 
transformation. 

Equations (3.8-3.11) connect the momentum variables via a Lorentz transformation that depends upon the inter- 
acting mass M\. The BT construction also uses a transformation that uses a non-interacting mass. To make this 
connection, we introduce constituent-particle momenta bi,b 2 ,b 3 . We use a different notation for these constituent 
momenta because they are not necessarily the same as Pi,P2,P3; for example, we will sec that b 3 ^ p 3 . We define 
variables analogous to Eqs. (3.8-3.10) for the non-interacting system 



(bi,b 3 ,P3)-»(B,b,k). 



(3.12) 
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The variables B and b are denned analogously to the definitions of P and p: 

B := bi2 + b 3 bi2 := bi + b 2 

b := A- 1 (B/M )6 3 ; b = b 3 + $_(b 3 , B, M )B, 



(3.13) 
(3.14) 



En 



Ev b; 



+ m 1 



Mn 



B 2 . 



(3.15) 



The additional variable k is obtained by (1) first boosting all three bi to the rest frame of the non-interacting three- 
body system, followed by a boost of particle 1 to the rest frame to the (12) pair. The vector k is the three-momentum 
of particle 1 after applying these two transformations. While the expression for k as a function of the momenta b^ 
can be written down explicitly, only the magnitude of k plays a non-trivial role in our calculations. It is related to 
the invariant mass of the non-interacting two and three-body systems by 



M fe := + k 2 M = yj M% + b 2 + ^/m 2 +b 2 . 

In the instant-form Bakamjian-Thomas construction, we identify 

B = P; b = p. 



(3.16) 



(3.17) 



The first of these identifications is a general property of an instant-form dynamics, while the second equation is a 
requirement of the Bakamjian-Thomas construction. 

These identifications and the above definitions imply relations between the observable momenta p 3 ,Pi2 and the 
constituent-particle momenta, b l7 b 2 , b 3 . It is because p and b are related to the momentum of particle 3 by Lorentz 
boosts involving different masses that p 3 ^ b 3 . 

The BT four-momentum eigenstate with particles 1 and 2 bound in the three-particle basis is 



(B, b, k|A; P', p% T = 5(P' - B)5(p' - b)0 A (k), 



(3.18) 



where <p\ (k) is the bound-state wave function. Note that b depends upon the free two-particle invariant mass Mi 20 
via Eq. (3.14-3.16), while p' depends upon the bound state mass eigenvalue via A Eq. (3.9-3.10). The association 
of p and b in the delta function of Eq. (3.18) is only consistent when P = 0. This is because in the instant form 
Bakamjian-Thomas construction the dynamical model is solved in the rest frame and the resulting mass eigenvalue 
determines how the system behaves under Lorentz boosts. The identification of p with b becomes inconsistent in 
frames where P / because these quantities are defined using different boosts. Since the initial and final states in 
a scalar density matrix element are in different frames, this inconsistency cannot be avoided and it results in the 
violation of cluster properties. 

We now make use of Eq. (3.18) to compute the BT counterpart to Ftp that was defined by Eq. (3.5): 



•Fbt ■= J dp' 12 ^Juj m (p' 3 ) B T(X, to; p' 12 , p 3 |j(0)|A, 

m t Pl2> P3) BT^J a -'ni(P3) 

:= / dp' 12 dk' dk^J uj m (p' 3 ) 



(3.19) 



d(P',p') 



<9(pi 2 ,P 3 ) 



d(p,p) 



<9(pi 2 ,p 3 ) 



d(b' 12 ,b' 3 ) 



d(P>,p>) 



<9(b i2 ,b 3 ) 



0(P,P) 



x 4>* x (k') (b' 12 , b 3 , k' \j (0) |bi2, b 3 , k)<h (k) ^JJj^j. 

The suffix k in the Jacobian indicates that the non-interacting two-body invariant mass, Mfe(Eq. 3.16), was used 
in the transformation of Eq. (3.14), while the suffix A indicates that the mass of the bound pair, M\ 1 was used to 
compute the variable change. The k-dependence of Jacobians with the subscript k is limited to the magnitude of k. 
Since the scalar density operator j(x) operates only in the space of particle 3, we have 



„(b' 3 )(b' 12 ,b 3 ,k'|j(0)|b 12 ,b 3 ,k)^ m (b 3 ) = mS(k' - k)<5(b' 12 - b 12 )G[(6 3 - b 3 ) 2 }. 
The integral over p' 12 can be converted to an integral over b' 12 . The Jacobian of the variable change is 



J dp' 12 = j dh' 12 



<9(P',p') 



9(b 12 /b' 3 ) 



<9(Pl2:P 3 ) 



5(P',p') 



(3.20) 



(3.21) 
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The final result is 



X"instant 
J- BT 



-I 



dk 



y/uj m (b' 3 ) 



d(Pl2,Ps) 



<?(P',p') 



0(P,p) 



9(P12,P3) 



d(p', P ') 



d(p, P ) 



\MV\' 



V^m(P3) 
V W ™( b 3) 



mG[(^-6 3 ) 2 ]- 



d(b' 12 ,b' 3 ) 

(3.22) 

The relevant observation is that this integral has a non-trivial dependence on the momentum, P12, of the bound 
state, in contrast to the P12 independence of the tensor product result. In this model the BT and tensor product 
representations of the three-particle system are related by one of Sokolov's unitary transformations, A, that preserve 
the three-body S'-matrix: 



A, to; p 12 , p 3 ) T p = A \ X > TO 5 P12, P3> 



BT- 



(3.23) 



The scale of the P12 dependence in Eq. (3.22) provides a measure of the size of the violation of cluster properties that 
results from ignoring A by replacing it with the identity. Since we know both states in Eq. ( 3.23) we do not have to 
calculate A explicitly to determine its impact. 

For this calculation, we vary the three-momentum transfer q and the spectator momentum pi 2 . The initial mo- 
mentum of particle 3 and the final momentum of the system are fixed in terms of these variables: 

• P3 = -gq; 

• P' = P + q. 

Since in the instant form, 



P = B: 



P' = B', 



(3.24) 



and the scalar density matrix element constrains b' 12 = bi 2 , we find that b' 3 = b 3 + q. However, the final momenta 
p' 12 and pg arc not constrained, and in general there are non-vanishing contributions to this matrix element for 
Ps P3 + q and p' 12 ^ pi 2 . 

In the nonrelativistic limit, where k, q and pi 2 are all small with respect to the relevant masses, the Jacobians are 
approximately unity and can be factored out of the integral, leaving a unit wave function normalization and a result 
identical to the TP case. The quantitative level of disagreement with the TP result is therefore linked to the extent 
to which the model goes beyond the nonrelativistic limit. 

For the other forms of dynamics we preform similar calculations. However, rather than using the natural variables 
for each given form of dynamics (four- velocity, light-front components of the four-momentum), we use the same 
variables employed in the instant-form case to facilitate a comparison of the size of the violations of cluster properties. 
These variables are also the quantities that are most readily measured using detectors. Calculations using front and 
point-form dynamics lead to results that have the same general form as Eq. (3.22) involving an integral over Jacobians 
and their inverses involving different sets of variables and the bound state wave functions. The Jacobians all cancel 
in the limit that two-body mass eigenvalue becomes the two-body invariant mass. 



D. BT Representation: Front Form 

Dirac's front-form dynamics is described in detail in Ref. [4]. In a front-form dynamics the generators of trans- 
formations that leave a plane tangent to the light cone invariant are free of interactions. We provide a summary 
here. 

Basis states in the front form are eigenstates of the light-front components of the momenta 

p = (p±,p+); P + =p°+p 3 ; P± = (p\ P 2 ). (3.25) 

These operators generate translations in a hyperplanc x~ := x° — x 3 = tangent to the light cone. In what follows 
we use a tilde to indicate the light-front components of four-vectors. 

In the front form the Lorentz transformations Eq. (3.9) and Eq. (3.14) used to define a particle momentum in the 
rest frame of the three-particle system in the instant form are replaced by boosts that leave the light front invariant. 

For the front-form Bakamjian-Thomas construction Eq. (3.17) is replaced by 

B = P; b = p. (3.26) 

As in the instant-form case the first equation is a general property of a light-front dynamics while the second equation 
is an additional requirement of the Bakamjian-Thomas construction. The second equation identifies the momentum 
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of particle 3 in the rest frame of the non-interacting three-particle system, p, with the momentum, b of particle 3 
in the rest frame of the (12) + 3 system. As in the instant-form case this identification has no consequences for the 
two-body S matrix, but will lead to a violation of cluster separability at the four-body level. With this modification 
the front-form BT expression for J^™* is similar to the instant-form result: 



x-front 
•> BT 




d(p'l 2 ,P 3 ) 
9(P',p') 



a(P,p) 

<9(P12,P3) 



d(P',p') 



d(bi2,b 3 ) 
0(P,p) 



i^(k)r 




6s) 2 ]. (3-27) 



This quantity has an unphysical dependence on bi2 that does not occur in the tensor product. 

The instant- and front-form calculations differ in both the constraints (3.26) vs (3.26) and the different boosts use 
to construct p, p',b and b'. For the purpose of the comparison with the instant-form calculations we make a final 
variable change replacing the light-front components of the momenta by the same variables that were used in the 
instant-form calculation. 

For this calculation, we vary q and pi 2 , with coordinate axes chosen such that q + = and pf 2 = 0. We define p_i_ 
and P' in terms of these quantities: 



P± 



• P' = P + q. 

Analogous to the discussion following Eq. (3.24), we have that 



P = B 



P' = B', 



(3.28) 



The scalar density matrix element constrains b' 12 = bi 2 , and therefore b 3 = b 3 + q. However, the final momenta p' 12 
and p 3 are not constrained, and in general the integral has non-zero contributions from p 3 ^ p 3 + q and p' 12 ^ pi 2 . 

The constraints pj_ = and P' = P + q ensure that the momentum transfer to the one-body system is q, as 

in the instant form example. 



E. BT Representation: Point Form 

Dirac's point-form dynamics are also described in detail in Ref. [4]. In this case the Lorentz group is non- interacting. 
We provide a summary here. 

Basis states in the point form are described by three-velocity vectors v. Momenta are obtained by multiplying the 
four-velocities by (interacting or non-interacting) masses. Thus, we seek to evaluate matrix elements of the scalar 
density operator j(x) between states of particle 3 and the (12) bound state with initial three-velocity V = P/M and 
final three-velocity V = P'/M'. 

In the point-form Bakamjian-Thomas construction, equations (3.17) and 3.26 are replaced by 



V = B/M ; b = p. 



(3.29) 



The first of these equations corresponds to total four-velocity conservation, which is a property of all point-form 
representations. The boosts in the point-form are the same as those used in the instant form, but they are param- 
eterized by the conserved four-velocities. The second equation, as in the instant and from form, is a requirement of 
the Bakamjian-Thomas construction. 

The derivations proceed in a fashion similar to the instant form. 

With these conventions, the point-form Tbt has the structure 



x-point 

BT -~ 



/^m(P 3 ) 

w m (b 3 ) 



<9(Pi 2 ,P 3 ) 



9(V',p') 



0(V,p) 



<9(P12,P3) 



0(V',p') 



d(bi2,b 3 ) 



0(V,p) 



i<Mk)r 



'wmCPjO 

w TO (b 3 ) 



mG[(b' 3 -b 3 n 



<9(b' 12 ,b 3 ) 

(3.30) 

Note in this case that the velocities V and V are the same for the kinematics using the interacting mass M or the 
free mass M . Since the boosts depend on V it follows that 



b 3 = p 3 ; b 3 = p 3 , 

and therefore that b 3 — b 3 = p 3 — p 3 = q. One consequence is that in general, M'V — MV ^ q. 
We constrain the momentum transfer to particle 3: 



(3.31) 
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• b' 3 = b 3 + q. 

These calculations share the property that b 3 — b 3 = q in the instant and point forms, and b 3 b 3 = q in the 
front form, i.e. the constituent-particle momenta are related to each other at the three-vector level by the momentum 
transfer to the overall three-body system. One could choose different sets of particle momenta in calculating BT 
matrix elements to compare to the TP result. For example, we could constrain the three calculations such that 
P' — P = q in some frame. Different choices explore different kinematic regions of the difference between A and the 
identity. Our investigations of other choices indicate numerical effects of comparable size. 



In this section we discuss results of calculations of the form factor Tbt of a scalar density for instant-, front- and 
point-form kinematic choices. Initially we consider scales that are relevant for systems of nucleons interacting with 
two-body interactions. Then we turn to examples more appropriate to hadronic models with subnucleon degrees of 
freedom. 

In all figures we display the figure of merit: 



which represents the relative error induced by ignoring the unitary transformation Eq. (3.23) that restores cluster 
separability. 



To model realistic conditions for nuclear physics, we use a bound-state spectator with deuteron properties con- 
structed from Malfliet-Tjon [22] potential IV, which contains both attractive and repulsive forces. Ref. [4] discusses 
how to make a phase-equivalent relativistic model that has the same wave functions as the non-relativistic model, and 
our numerical results follow that approach. The invariant form factor G(q 2 ) is taken as a dipole form factor. 

We first consider the figure of merit as a function of the momentum transfer and momentum of the bound state in 
each of Dirac's forms of dynamics and for the bound state momentum perpendicular and parallel to the momentum 
transfer. The results are shown in figures 2-7. In all cases the expected results are given by the flat planes. 

The fractional deviation of this BT model calculation from the TP result that satisfies cluster separability is very 
small for all three forms of dynamics, typically of order 10 -3 or smaller at the highest values of the three-momenta q 
and pi 2 . 

instant form (BT-TP)/TP vs. q, P 12 ; q parallel to P 12 instant form (BT-TP)/TP vs. q, P 12 ; q perp to P 12 



IV. RESULTS AND DISCUSSION 



(J~BT — J~tp) 

Ftp 



(4.1) 



A. Malfliet-Tjon Deuteron Wave Function 



1.5 x 10' 3 

1.0 X lO 3 
5.0 X 10~ 4 
0.0 X 10° 




FIG. 2: (Color online) Model differences for instant-form 
BT calculation, q||pi2- 



FIG. 3: (Color online) Model differences for instant-form 
BT calculation, q _L pi2. 
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front form (BT-TP)/TP vs. q, P 12 ; q parallel to P 12 



front form (BT-TP)/TP vs. q, P 12 ; q perp to P 12 
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FIG. 4: (Color online) Model differences for front-form 
BT calculation, q||pi2- 



FIG. 5: (Color online) Model differences for front-form 
BT calculation, q _L pi2. 



point form - (BT-TP)/TP vs. q, P 12 ; q parallel to P 12 



point form - (BT-TP)/TP vs. q, P 12 ; q perp to P 12 




P12 [fm" 1 ] 




FIG. 6: (Color online) Model differences for point-form 
BT calculation A, q||pi2- 



FIG. 7: (Color online) Model differences for point-form 
BT calculation A, q _L pi2. 



B. Binding Energy Variation 

The calculations above assumed a bound state with a wave function having a typical dependence on the relative 
momentum of the constituent nucleons. The next set of curves illustrates the hgure of merit for fixed values of q and 
P12 as we vary the binding energy and momentum scale of the wave function. The variations that we consider are 
still scales that are relevant to nuclear systems. Figures 8 and 9 show the results of calculations for varying binding 
energy with a Malfliet-Tjon wave function. 

As with the earlier cases that employed the bound-state binding energy, the fractional deviation of the BT results 
from the TP benchmark is quite small, of order 10 -3 or less. 



Wave Function Scale Variation 



We also examined sensitivity to the momentum scale, fcrj, of the wave function by replacing the Malfliet-Tjon wave 
function with a Gaussian form: 



i 



,-(fc/feo) 2 



N 



(4.2) 



Figures 10 and 11 show the results for the spectator momentum perpendicular and parallel to the momentum transfer 
in all three forms of dynamics. 

These results mirror those discussed above: the fractional deviation of the BT results from the TP benchmark is 
quite small, of order 10 -3 or less. 
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all forms - (BT-TP)/TP vs. BE, q = 10, P 12 parallel = 5 



all forms - (BT-TP)/TP vs. BE, q = 10, P 12 perp=5 
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FIG. 8: (Color online) Model differences for BT calcula- 
tions as a function of two-body spectator binding energy, 
q=10fm- 1 ||p 12 . 



FIG. 9: (Color online) Model differences for BT calcula- 
tions as a function of two-body spectator binding energy, 
q = 10 fm -1 _L p i2 . 



all forms - (BT-TP)/TP vs. kO, q = 10, P 12 parallel=5 
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FIG. 10: (Color online) Model differences for BT calcu- 
lations as a function of two-body Gaussian wave function 
scale, q = 10 fm _1 |]pi2. 



FIG. 11: (Color online) Model differences for BT calcu- 
lations as a function of two-body Gaussian wave function 
scale, q — 10 fm -1 _L pi2. 



D. Implications for Nuclear Theory 

Formally dynamical relativistic models should be Poincare invariant and satisfy cluster properties; in addition 
operators associated with short ranged physical phenomena should have well defined cluster expansions. 

While our model is considerably simpler than a realistic model of a three-nucleon system interacting with an 
electron, the difference between the Bakamjian-Thomas formulation of the model and the model that clusters properly 
is due entirely to the same unitary transformations discussed in our simplified model. This difference vanishes in the 
limit that these transformations become the identity. The analysis in this paper established that these operators, 
which implicitly depend on the interactions, are close to the identity for interactions with typical nuclear physics 
scales for binding energy and Fermi momentum. This provides a strong justification for applying Bakamjian-Thomas 
construction of dynamical representations Poincare group to system of more than three nucleons. 

Our results indicate that Bakamjian-Thomas models, which explicitly satisfy the requirements of Poincare invari- 
ance, can be utilized in typical nuclear physics problems with minimal quantitative error due to the lack of cluster 
separability using any of Dirac's front-form dynamics. 



E. Implications for Hadron Models 



The final set of figures show the results of calculations with scales that are more appropriate models of hadrons 
based on sub-nucleonic degrees of freedom. 

We note here that confinement precludes separating arbitrary subsystems by large distance scales, so the requirement 
of cluster separability is irrelevant for models of individual hadrons. The issue may be relevant, however, for systems 
of hadrons described by sub-nucleonic degrees of freedom. 
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In these cases, in order to understand the relevant scales, we replace the particle 3 mass in the above calculations 
by a "constituent quark" mass of 220 MeV, and consider two-body ("diquark") masses ranging from 200 to 600 MeV 
for Gaussian wave functions with a 1 fm" 1 scale in Figs. 12 and 13, and Gaussian wave functions with scale ranging 
from 0.5 fm -1 to 10 fm -1 , with a diquark mass of 600 MeV in Figs. 14 and 15. 

The results are similar to those appropriate to nuclear physics discussed above, except that the scale of deviation 
from the TP benchmark is on the scale of 10% for the instant, front and point-form As. This not surprising given that 
the mass/momentum scale variation for these calculations is much higher than for typical cases in nuclear physics with 
nucleons. While these corrections may be relevant, multi-hadron models based on sub-nucleon degrees of freedom 
have not been developed to a precision where these effects could be identified. 



all forms - |(BT-TP)/TP| vs. md, q = 10, P 12 perp=5 



all forms - |(BT-TP)/TP| vs. md, q = 10, P 12 perp=5 



b. 



1 x IO" 3 

1 X 10" 

l x icr 5 
l x io 6 

1 x 10" 7 





















instant perp \ 




front perp - - 




point perp 



300 350 400 450 500 

diquark mass [MeV] 



550 600 



1 x 10" 1 

1 X IO 2 

t. i x 10 3 

Q_ 

fe 1* 10 4 
1 X IO 5 

1 x IO" 6 







instant perp - 




front perp - - 




point B perp - 



200 250 300 350 400 450 500 550 600 

diquark mass [MeV] 



FIG. 12: (Color online) Model differences for BT calcu- 
lations as a function of two-body spectator mass, q 
10 fin-^lpia. 



FIG. 13: (Color online) Model differences for BT calcu- 
lations as a function of two-body spectator mass, q = 
10 fm" 1 _L pi 2 . 
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FIG. 14: (Color online) Model differences for BT cal- 
culations as a function of Gaussian wave function scale, 
q = 10 fm _1 ||pi2. 



FIG. 15: (Color online) Model differences for BT cal- 
culations as a function of Gaussian wave function scale, 
q = 10 fm -1 _L pi2. 



V. SUMMARY 



Bakamjian- Thomas formulations, which explicitly satisfy the requirements of Poincare invariance, do not satisfy 
cluster separability above the three-particle level, i.e. in systems that involve three-body subsystems whose total 
momentum must vary. The cluster properties can be restored via a hierarchy of unitary transformations. These 
transformations depend upon the full solution of a three-body problem, and are difficult to implement in practice. 

Rather than attempt to calculate directly the size of these unitary transformations (e.g. the difference of matrix 
elements from those of the unit operator) , we have developed a simple model in which the exact result consistent with 
cluster separability is known, and then compare to it the results of Bakamjian-Thomas calculations. The quantities 
that we investigate are sensitive to the part of these unitary transformations that restore cluster properties. 
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We conclude from these model studies that Bakamjian-Thomas models, which explicitly satisfy the requirements 
of Poincare invariance, can be utilized in typical nuclear physics problems with minimal quantitative error due to the 
lack of cluster separability. 

We also examined models utilizing mass/momentum scales appropriate for quark models. Confinement precludes 
separating arbitrary subsystems by large distance scales, so the general requirement of cluster separability is irrelevant 
for models of individual hadrons. The issue may be relevant, however, for systems of hadrons described by subnuclconic 
degrees of freedom. In such cases, the deviations from the model benchmark are larger (about 10%) than those for 
typical nuclear physics calculations with nucleons, though they are still small compared to model uncertainties. 

This work was supported in part by the U. S. Department of Energy, Office of Nuclear Physics, under contract No. 
DE-FG02-86ER40286. 
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